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NON-LINEAR INTERACTION OF A 
DETON ATION/VORTICITY WAVE 1 

D. G. Lasseigne & T. L. Jackson 
Department of Mathematics and Statistics 
Old Dominion University 
Norfolk, Virginia 23529 

M. Y. Hussaini 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, Virginia 23665 


ABSTRACT 

The interaction of an oblique, overdriven detonation wave with a vorticity disturbance is 
investigated by a direct two-dimensional numerical simulation using a multi-domain, finite- 
difference solution of the compressible Euler equations. The results are compared to those 
of linear theory, which predict that the effect of exothermicity on the interaction is relatively 
small except possibly near a critical angle where linear theory no longer holds. It is found 
that the steady-state computational results whenever obtained in this study agree with 
the results of linear theory. However, for cases with incident angle near the critical angle, 
moderate disturbance amplitudes, and/or sudden transient encounter with a disturbance, 
the effects of exothermicity is more pronounced than predicted by linear theory. Finally, it 
is found that linear theory correctly determines the critical angle. 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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I INTRODUCTION 

The passage of a weak shear disturbance through a reactive shock wave, or deto- 
nation, was examined by Jackson, Kapila, and Hussaini^. Supersonic engines based on 
oblique, overdriven, reacting shock waves have been proposed as a possible alternative 
to the SCRAMJET for high-speed propulsion 2 . It is still a mattter of exploration as to 
whether or not such waves can be stabilized. Of particular interest in this investigation was 
the effect of heat release on the refraction and amplification of the vorticity disturbance 
and the simultaneous generation of acoustic and entropy signals behind the overdriven 
detonation. The detonation was assumed to be at an angle to the base flow, and the nor- 
mal Mach number of the gas ahead of the detonation front was taken to be greater than 
the Mach number of a Chapman-Jouget wave. The. vorticity disturbance was assumed to 
be a small amplitude, planar, shear wave with wave vector parallel to the base flow (i.e. 
transverse disturbances). There exists a critical angle, dependent upon the exothermicity 
of reaction and the overdrive, such that the relative velocity of the base flow behind the 
front is subsonic for 9 < 9 C and supersonic for 9 > 8 e . In the former case, the amplitude 
of the generated acoustic disturbance is exponentially decaying behind the detonation; 
in the latter case, the amplitude is constant. The critical angle approaches zero as the 
Chapman-Jouget limit is approached. It was found that the vorticity was significantly 
amplified by the exothermicity and that the generated acoustic response is most affected 
by exothermicity near the critical angle. Furthermore, the manner in which the shape and 
structure of the detonation are altered by the disturbance was also investigated. 

The analysis was accomplished by considering the detonation wave, which consists 
of a lead shock, an induction zone and a fire zone, as a discontinuity on the length scale 
of the disturbance. The discontinuity separates an unburnt mixture of reactants from a 
burnt mixture of reaction products. After superimposing a vorticity disturbance in the 
form of a planar shear wave, the generalized Rankine-Hugoniot conditions were linearized 
about the base flow providing the conditions to determine the amplitudes and angles of 
the transmitted vorticity wave and the generated acoustic and entropy waves. The effect 
of the disturbance on the structure of the detonation was found by considering the limit 
of large activation energy, E. The thickness of the induction zone (the region between the 
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lead shock and the fire zone) is a measure of the thickness of the detonation. The lead 
shock and the fire zone were treated as discontinuities on a length scale comparable to the 
induction zone thickness. The equations governing the perturbations within the induction 
zone due to disturbances ahead of the lead shock were derived, and the fire zone position 
was determined as the point at which the the solutions become singular. It was found that 
disturbance amplitudes 0(E~ 1 ) have an 0(1) effect on the fire zone position. 

The aim of the present investigation is to examine the predictions of linear theory 
and to determine by means of a direct two-dimensional numerical simulation the regions of 
validity for linear theory. For this investigation, the detonation is treated as a discontinuity, 
and the generalized Rankine-Hugoniot relations are used to provide appropriate jump 
conditions across the detonation. The effect of the disturbances on the internal structure 
of the detonation is not considered. Of particular interest is the behavior of the solution 
for angles near critical where linear theory, if it applies, predicts that exothermicity can 
have a significant effect on the generation of acoustic signals. 

II NUMERICAL SIMULATION 

To study the interaction of a detonation wave with a vorticity disturbance, a numerical 
approach similar to the one utilized by Zang, Hussaini, and Bushnell 3,4 is considered. 
The major differences are that the present routine takes advantage of the periodic nature 
of the solution in the transverse direction and also uses a multi-domain scheme in the 
normal direction. Rather than using coordinates in which the position of the detonation 
is stationary, it is assumed that at time t = 0 an infinite, planar detonation wave starts at 
x = 0 and propagates into the unburnt mixture. The position of the detonation x,(y,t ) is 
calculated by using the improved shock-fitting approach of Kopriva, Zang and Hussaini 5 , 
generalized to include exothermicity. The initial conditions are chosen such that in the 
absence of disturbances the detonation wave will propagate to the right with a Mach 
number (relative to the gas in front of the detonation wave) greater than the Chapman- 
Jouget Mach number. The base flow in the region ahead of the detonation has only a 
vertical component so that the net result of the velocity of the detonation front in the 
positive x- direction and the vertical component produce a base flow at an angle 0 in a 
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coordinate system attached to the detonation front. Superimposed on the uniform base 
flow is a vorticity disturbance which propagates at the same angle 6 to the x-axis (see Figure 
1). The disturbance is assumed to have unit wavelength, 27r/|fc| = 1 where k is the wave 
vector. The flow is therefore periodic in the y-direction with period 27 t /k y = sec 6. The 
initial flow behind the detonation front has a vertical component such that the tangential 
velocity across the front is continuous. It should be noted that the flow ahead of the front is 
prescribed and is used to impose the appropriate jump conditions across the discontinuity. 

The physical domain in which the fluid motion is computed is given by 

x L < x < x,(y,t), 0 < y < sec(0), t > 0, (1) 


where x, is the shock position and x, (y, 0) = 0. The left boundary x L is some suitably 
chosen negative number (usually minus one). This domain takes advantage of the periodic 
nature of the solution in the y-direction. In the x-direction, a multi-domain approach is 
used which allows for a greater number of grid points near the detonation front where 
greater accuracy is needed and fewer grid points away from the front where the solution 
is smoother. The interfaces between the domains are denoted by z*(y,t ) where the left 
boundary is xo(y,t) = x L , the detonation front is x n (y,t ) = x ,{y,t), and the planar 
interfaces between the domains are Xi(y,t) = x;(t),fori = l,n-l. Since the behavior of the 
solution near the critical angle is of interest, it is necessary to continue the calculations until 
the detonation wave penetrates a significant distance into the unburnt mixture. Because 
of this expanding domain, a greater number of grid points in the x-direction is needed 
at the end of the calculation than is needed in the beginning. The resolution in each 
subdomain is monitored, and the number of grid points is increased when necessary. A 
linear interpolation is used to transfer the solution to the finer grid. The numerical routine 
also allows for ramping of the amplitude of the imposed disturbance, thereby decreasing 
in size the portion of the transient caused by the sudden encounter of the undisturbed 
detonation with the vorticity wave. The approach to steady state was smoother, allowing 


an accurate comparison to linear theory. 

The following change of variables is made to computational coordinates 

x _ — x - ; _ i, „ 

•£«(j/i^) i(y> 0 


(2.a) 
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Y = y cos(0) 
T = t 


( 2 . 6 ) 

(2.o) 


The computational domains are therefore 


0 < X t < 1, 0 < Y < 1, T > 0. 


(3) 


The fluid motion is governed by the two-dimensional Euler equations. In terms of the 
computational coordinates these are 


where 


Qt 4- BQ Xi + CQy =0 i = 1, n 


0 = [P^.UjSf 


B = 


Ui iX im ~fXi. 


0 1 


c 2 X ix /i Ui 0 0 

0 Ui 0 


c 2 x iy h 


c = 


0 0 0 Ui 

V iY x 7 r y 0 

c 2 y, /7 v o o 


c 2 y„/ 7 o 


o 


o o o y J 

Ui = Xi t -f- llXix + VXiy 


v = y t + u y« + vY y 


(4) 


(5a) 


(56) 


(5c) 


(5 d) 
M 


Here, P is the natural logarithm of pressure, c is the local sound speed, and 5 is the 
entropy divided by specific heat at constant volume. The velocities u and v, in the x and 
y directions respectively, are scaled by the sound speed ahead of the shock. The ratio of 
specific heats 7 is taken to be 1.4 for all calculations. 

As in Zang, Hussaini and Bushnell 3,4 , the equations are discretized in each subdomain 
using the finite difference method of MacCormack 6 , which is a variant of the Lax-Wendroff 
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method. Conditions at the right boundary, at the interface between the domains, and at 
the left boundary need to be imposed to determine the system completely . At the right 
boundary, the improved shock-fitting routine of Kopriva, Zang and Hussaini 7 has been 
generalized to include exothermicity and time dependence of the flow ahead of the shock. 
Details are given in the Appendix. At the interface between domains, a routine similar to 
Kopriva 8 is used. Since the interfaces are perpendicular to the A^-axis, the derivatives in Y 
can be calculated by using the values of Q along the interface. The values of the derivatives 
in X , however, must be taken from the left or right of the interface or a combination of 
both. The method chosen here is 


o; + i (b + |flr)<R + \(b - |b|-)q; + cq' = o (6) 

where \B\* is an approximation to |i?| = P|A|P -1 , A is a diagonal matrix with the eigen- 
values of B on the diagonal and P is a matrix of the eigenvectors of B. The superscript I 
refers to values calculated on the interface and the superscripts R and L refer to the finite 


difference approximations to the X derivative in the subdomain to the right and to the 
left of the interface respectively. The approximation used for these calculations is 
|" \\U + c<j>\ + \\U - c<f>\ 0 \\U -I- c<j>\ — \\U — C(f)\ 0 

0 |c/j o o 

\\U + c<t>\ - \\u - c<t>\ 0 \\XJ + c4>\ 4- \\U — ccj)\ 0 

o o o \u\-> 


\B \ * = 


( 7 ) 


where <p = (X*+Xf) 1/2 . The conditions imposed at the left boundary naturally depend on 
the inflow at the boundary being supersonic or subsonic. If the inflow is supersonic, then 
all four components of Q can be prescribed; however, if the inflow is subsonic, only three 
of the four components can be prescribed at the left boundary. For the parameters of this 
study, the inflow is always supersonic. It should be mentioned that the normal component 
of the flow behind the detonation is supersonic in the laboratory frame; however, the 
normal component of the flow behind the detonation is subsonic in a frame of reference 
moving with the detonation. 
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Ill RESULTS AND COMPARISON WITH LINEAR THEORY 


Linear Theory. 

Consider a coordinate system attached to the detonation front such that the speed of 
the detonation is zero. The flow ahead of the front is assumed to be at an angle 0 to the 
front, the normal component of the velocity has Mach number Mq, and the vertical com- 
ponent of the velocity is Wo (see Figure 1). The generalized Rankine-Hugoniot conditions 
for a detonation with heat release parameter a provides that the normal component of the 
flow behind the front has Mach number (relative to the upstream speed of sound) 


mo = 


1 -I- 7 M$ 


(Ml ~ 1 ) ; 

~~Mf~ 


- 2 ( 7+1 )a 


-i 1/2 


(7 + l)Mo (7 + 1) 

and requires that the tangential components be continuous 

wo = Wq, 

Therefore, the angle of the base flow behind the detonation is given by 

M 0 

tan <p = tan 0, 

mo 


( 8 ) 


( 9 ) 


( 10 ) 


Superimposed on the upstream base flow is a sinusoidal disturbance in the velocity given 
by 

U\ = e&o cos(fc' • £) (11) 

where Vq = ( Mq,Wq ) and Vq • k' = 0. When this constant pattern of vorticity is con- 
vected through the detonation front, the vorticity wave is refracted and amplified, and an 
acoustic and entropy wave are generated. Linear analysis provides analytic expressions for 
the amplitudes and angles of all three disturbances downstream. The result is that the 
following planar waves are superimposed on the base flow and the pressure downstream 
( x < 0 ): 


ui = 


U 1 

Vl 


= e|tfo|S 


cos(fc • x + 6 , ) 
0 


+ e\U 0 \P(x) 


cos (k p • £ + 6 p ) 
Psin (k p ■ x + <5 p )J 


( 12 ) 


PI = — e7|Z7o|Mo sec<j>P(x) cos (Jc p ■ x + 6 p ) 


(13) 
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where 


( 5° y > 1 
S={ 


(14) 


l So A 4 ^ 1) 

r p° n > i 

P = { (15) 

l Poexp(Aa;) n < 1, 

and ^ is the local Mach number of the relative base flow behind the detonation front. The 
constants 5°, So, P°, Po, <5,, S, and A are given in Jackson, Kapila, and Hussaini 1 . The 
vector k is perpendicular to the downstream base flow and the vector k p is perpendicular to 
the direction of acoustic propagation which makes an angle <j>' with the x-axis (see Figure 
2). This angle is given by 




4> — cot -1 (/?) 


n > 1 


(16) 


l —tan 1 [(^//?„) 2 cos 2 </>tan<^] (j, < 1 

with (3 = (|l — //^|)1/ 2 and /? 2 = 1 — /i 2 cos 2 A similar disturbance is also superimposed 
on the temperature and density due to the generated entropy wave. 


Nonlinear Calculations. 

The comparison of the two-dimensional calculations to the linear theory consists of 
comparing the amplitudes of the vorticity and acoustic waves behind the detonations with 
the predictions. To achieve this comparison, a least squares fit to the functional forms 

u y - v a = Sicos(fc'y) + 52 sin(fc' y) (17) 

P = Po + Pi cos(fc^y) + P 2 sin(k' y y) (18) 

is performed for the computed vorticity and acoustic fields for each value of x. The 
amplitudes of the above vorticity and acoustic disturbances 5 = (S 2 + 5 2 ) 1 ^ 2 and P = 
(P x 2 + Pi) 1 / 2 are divided by e\Uo\k' y and ei\Uo\Mosec<j> respectively and graphed as a 
function of x for various times t. In the calculations presented here, the front is initially 
at position x = 0, and behind the front there is no disturbance (i.e. P = 0 and 5 = 0). 
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As time progresses, the front moves downstream into a region with a non-zero vorticity 
disturbance. Therefore, the numerical solutions include transient behavior. The relaxation 
to steady state depends on both the acoustic and the vorticity responses as well as the 
angle of incidence. If 6 is above the critical angle, then the predicted acoustic response 
is constant; therefore, steady state is indicated by a broad flat response behind the front. 
If 9 is below the critical angle, an exponentially decaying acoustic response is predicted; 
therefore, steady state is indicated when the value of the acoustic response at the front 
remains essentially constant as the front propagates a sufficient distance. Since the vorticity 
response is constant in both cases, steady state is indicated when a broad flat vorticity 
response behind the shock has been achieved. 

For many of the calculations involving either large disturbance amplitudes or angles 
of incidence near critical, the solution did not asymptote to a steady state even after the 
detonation front had penetrated a significant distance into the unburnt mixture. The 
nature of this transient behavior is important since real turbulence consists of sudden 
non-steady phenomena, not steady plane waves. In general, the transient depends on 
the suddenness in which the front meets with the full disturbance, the amplitude of the 
disturbance, and the angle of incidence, TTo resolve the portion of the transient caused by 
the sudden encounter of the front with a non-zero disturbance, the calculations allow for 
a slow or fast ramping of the amplitude of the imposed disturbance; the amplitude as a 
function of time is given by 


e(t) = 


e 0 (M 0 t/R) 2 (3 - 2 M 0 t/R) 

eo 


t < R/Mo 
t > R/Mq. 


(19) 


Therefore, when the position of the front is x, = R, the flow ahead of the front will have 
its maximum disturbance amplitude. 

The value of the exothermicity parameter a is chosen to be equal to two for all 
calculations. This value is chosen since linear theory predicts that the value of 5 will be 
significantly different from its value in the non-reacting case for almost all values of angle 
9 and the value of P will be significantly different for values of the angle near critical. 
The strength of the detonation is chosen so that the normal Mach number of the flow is 
1.5 times the Chapman- Jouget number. The critical angle for these parameter values is 
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9 e = 24.89° 

The pressure and vorticity responses to a 1% disturbance at an angle of 9 = 40 for 
ramping parameter R = 6.0 and R = 0.5 are shown in Figures 3a and 3b, respectively. For 
these calculations, 32 grid points in the ^-direction are used. For slow ramping, the gas 
behind the front responds to each small increment in the amplitude of the disturbance, 
and a smooth approach to steady state is observed. For fast ramping, R = .5, the acoustic 
response shows significant overshoot of the predicted linear value. The vorticity response is 
somewhat smoother. The influence of heat release on the transient is considered by letting 
a = 0, R = .5 (Figure 3c) and prescribing the same normal Mach number; the overshoot 
and fluctuations of the acoustic response are considerably less than for the reacting shock, 
but the vorticity response shows essentially the same behavior as in the detonation. 

As seen from Figure 3a, when the ramping is slow and the amplitude of the disturbance 
is small, the predictions of linear theory agree with the calculated acoustic response but 
there is a discrepancy in the vorticity response of about 2%. As mentioned in Zang, 
Hussaini & Bushnell 3 ' 4 , since the vorticity response involves computation of derivatives in 
both the x and y directions, the vorticity response calculations are possibly less accurate 
than the pressure response calculations. When only 16 grid points are used in the y 
direction, the discrepancy in the vorticity response increases to 4-5% and the pressure 
response still agrees with the linear predictions. It should be noted that the calculated 
vorticity response is always less than the predicted values. Another important parameter 
which influences the transient behavior is the amplitude of the disturbance. For the 10% 
disturbance at 9 = 40° (Figure 3d), the acoustic response takes a longer time to asymptote 
to a steady state. 

The comparison between linear theory and numerical simulations for other angles is 
seen in Figure 4. The circles are the calculated responses to a 1% disturbance and the 
diamonds are the calculated responses to a 10% disturbance. As previously found for 
the interaction of a non-reacting shock with a vorticity disturbance, the linear predictions 
and the calculated acoustic responses agree for disturbance amplitudes up to 10% when 
the angle of the incoming flow is not near the critical angle. The calculated vorticity 
responses are consistently about 2% below the predicted values when 32 grid points are 
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used. Although many calculations were made for small amplitude disturbances (1% or 
less) and angles within 5° of critical, the results are not given in Figure 4 since a steady 
state value could not be reliably determined. The results of these calculations are discussed 
below. 

For disturbance amplitudes larger than 10%, the solution did not reach a steady 
state. A typical run is shown in Figure 5 for a 30% disturbance at 9 = 40° with ramping 
parameter R = 3. Linear theory is a good predictor of the scale for the overall response, 
but the distortions produced in the detonation front prevent the solution from reaching 
steady state; in fact, the distortions prevented the calculations from proceeding further 
since the time step based on the CFL number became exceedingly small. Furthermore, 
since the front is no longer planar, the comparison between the least squares fit to the 
forms (17) and the linear case no longer holds. The same behavior was observed for 10% 
disturbances at angles within 5° of critical. 

Linear theory predicts that exothermicity has its greatest effect on both the vorticity 
and pressure responses near the critical angle; therefore, it is important to resolve the 
behavior of the solution for angles near critical. As the angle of incidence approaches the 
critical angle, it is observed that the relaxation of the transients takes longer and longer. 
These transients are not related to the ramping of the disturbance amplitude as was shown 
for 9 = 40°. The response to a 1% disturbance with 9 = 30° is shown in Figure 6 (16 
grid points in the y-direction and ramping parameter R = 6). The transient nature of the 
solution provides for considerable overshoot of the predicted acoustic response with very 
little relaxation over the duration of the run. To continue the calculations past x , = 10 
would require adding additional grid points in the x direction to keep the same resolution 
making each time step considerably more expensive. In addition to the large pressure 
fluctuations, the vorticity response also displays some overshoot which is not observed 
for slowly ramped disturbances at angles more than 5° from critical. This is particularly 
surprising since calculations at other angles indicate that the vorticity response remains at 
5% below the predicted value and it is quite smooth for such slowly ramped disturbances. 
The same parameter values were used but with fewer grid points in the ^-direction allowing 
integration to continue until x, = 24. Although the resolution is less, a relaxation of the 
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transients is observed; the acoustic response oscillates about a value accurately predicted 
by linear theory, and the vorticity response oscillates about a value which is approximately 
5% below the value predicted by linear theory; there appears to be a slow decrease in the 
amplitude of these fluctuations. These results are consistent with steady states at larger 
angles of attack found using 16 grid points. A similar trend in the vorticity response for 
angles of attack slightly less than critical develops. The vorticity response displays a slight 
overshoot of the value predicted by linear theory which is not present for angles away 
from critical, and returns to a small oscillation about a value below the linearly predicted 
value. The overshoot of the predicted pressure and vorticity responses is observed for 
angles of attack between two and five degrees of the critical angle; however, relaxation of 
the transients is significantly slower as the angle of incidence nears critical. For 9 — 24.9 
(Figure 7), it can be seen that the acoustic response does not reach the value predicted by 
linear theory but is still increasing when the calculation is stopped. The vorticity response 
is also still increasing, having reached only half of the value predicted by linear theory. 
This calculation was terminated due to the increasing expense of maintaining an accurate 
numerical resolution; however, the calculations were repeated using a very coarse grid, 
and it was found that even though the pressure and vorticity responses continue the slow 
increase, the pressure has achieved only 40% of the linear theory value and the vorticity 
only 60% of the linear theory value when x, =50. These results are only qualitative at 
best but seem to indicate that the lengthening of the transient response for angles within 
2° of critical counteracts the transient overshoot of the pressure and vorticity responses 
found for other angles within 5° of critical. 

IV CONCLUSIONS 

Nonlinear calculations of the response of an initially plane detonation wave to a vor- 
ticity disturbance show that the results of steady-state linear theory are useful in providing 
an overall scale of the response. In cases in which the angle of incidence is near critical, 
disturbance amplitudes are moderate, and/or there is a sudden encounter with a distur- 
bance, the calculated responses display a transient overshoot of the linear prediction. It is 
found that exothermicity increases the overshoot in addition to increasing the value of the 
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predicted linear response. Also, the significant departures of the predicted responses of 
the reacting shock from the non-reacting shock near the critical angles appear to be real, 
and the calculated responses show transient overshoot of the predicted values for angles 
between two and five degrees of critical. Closer to the critical angle, the lengthening of the 
relaxation time for the transient produces a competing effect. Also, it is found that criti- 
cal angle of linear theory is an accurate predictor of the transition in the behavior of the 
acoustic response. Previously 8 '^, it was reported that the change from a constant pressure 
response to an exponentially decaying pressure response occurred at an angle significantly 
different than predicted by linear theory. For the situation presented here, it is determined 
that the long relaxation of the transients for angles slightly above critical makes the pres- 
sure response appear to be exponentially decaying; when in fact, if the calculations are 
of a sufficient duration, the pressure response eventually changes and eventually oscillates 
about a a constant value which is consistent with the value predicted by linear theory. 

Since exothermicity is seen to increase both the vorticity and acoustic responses for 
the detonation wave/vorticity wave interaction, it would be appropriate to study the in- 
teraction of a detonation wave with a fully developed vortex. Such a numerical study was 
presented in Meadows, Kumar and Hussaini 8 for the non-reacting shock using a shock 
capturing scheme which is more appropriate than the shock fitting scheme for the study of 
shock interaction with large disturbances. Of particular interest is the ability to capture 
secondary shocks which agree with experimentally observed features of the flow. This is 
not possible using shock-fitting methods. 
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APPENDIX: 

The improved shock fitting scheme utilized by Kopriva, Zang, and Hussaini 5 is gener- 
alized to allow for exothermicity of the reacting shock and to allow for a time dependent 
flow ahead of the front. The generalized Rankine-Hugoniot conditions for a stationary 
reacting shock can be simply written as 


P2 = Pi + g{M) 


m = f(M ) 

M(1 + a + ^^(M 2 — m 2 )) 

m 


g(M) - In 


( 20 ) 

( 21 ) 

( 22 ) 


where the subscript 1 refers to known quantities ahead of the front and subscript 2 refers 
to the corresponding quantity behind the front; /(M) is given by (8). M is the normal 
Mach number of the flow ahead of the front, and m is the normal velocity behind the front 
divided by the sound speed ahead of the front. If all other velocities are nondimensionalized 
with respect to the speed of sound for the mean flow ahead of the front, then the prescribed 
flow has a local speed of sound c\ at the front which differs from unity and depends upon 
both y and t. Denoting the front position by f) (y, t) , the speed of the front traveling in 
the positive x direction is (y , £) = (Uj , 0). The velocities of the gas in the frame of 
reference where the front moves = (tx, , i \ ) are related to M and m by 


c\M = U f N, (23) 

c\m = Uj N m — Q 2 • $ (24) 


where i'? = ( N B , N y ) is the normal to the shock front pointing in the direction of the domain 
subscripted by 1. The key to the shock fitting routine is the compatibility equation, 

P 2i + • &>t = —3 • VP - P. - ctf ■ VP - R P (25) 

c c 

where 

R — (uu, -(- vu v , uv t + vv y ) Rp = 7 u, -I- 7 u„ (26) 

derived from the Euler equations. By differentiating the Rankine-Hugoniot relations and 
equations (23)and (24), an equation for the acceleration of the detonation front can be 
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found and is given by 


U }i (y,t) = {-cicPi, + c\cC + (cG - 'yF)(Qi i -I? + Qi • 

-(cG - 7-F + JV, t + C17Q2 • (27) 

+ci , [(cG - jF)M + ci7m]} {(cG - 7F -f 7Ci)iV s }~ 1 

where G = g'(M), F = f'(M), and C is the right hand side of equation (25). During 
the calculations, the quantity C is evaluated using the solution of the Euler equations 
from the previous time step, and all other quantities are calculated using the appropriate 
jump conditions. The front position and velocity are updated for each time step using the 
two-step Lax-WendrofF scheme. 
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Figure 1. Schematic of model problem in the physical domain. Velocity of unburnt 
mixture is Wq + IVj. 
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Figure 2. Schematic of model problem for linear theory. Acoustic disturbances in 
the burnt gas propagate at angle 4>' to the x-axis. Vorticity and Entropy disturbances are 
convected with the mean flow at angle <f>. 
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Figure 3d. Response to a 10% disturbance with slow ramping. R = 6, a = 2, and 
40®. The curves correspond to t = .237, .612, .704, .933, 1.208, 1.342* 
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Figure 4. Pressure and vorticity responses as a function of angle. Circles represent 
1% disturbances and diamonds 10% disturbances. Solid line is linear theory for a = 2. 
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Figure 7. Transient response to a 1% disturbance for 9 = 24.9°. R 
The curves correspond to t = ,562, .917, 1.335, 1.639, 1.999. 
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